Correlated motion of two atoms trapped in a single mode cavity field 
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We study the motion of two atoms trapped at distant positions in the field of a driven standing 
wave high Q optical resonator. Even without any direct atom-atom interaction the atoms are coupled 
through their position dependent influence on the intracavity field. For sufficiently good trapping and 
low cavity losses the atomic motion becomes significantly correlated and the two particles oscillate 
in their wells preferentially with a 90° relative phase shift. The onset of correlations seriously limits 
cavity cooling efficiency, raising the achievable temperature to the Doppler limit. The physical 
origin of the correlation can be traced back to a cavity mediated cross-friction, i.e. a friction force 
on one particle depending on the velocity of the second particle. Choosing appropriate operating 
conditions allows for engineering these long range correlations. In addition this cross-friction effect 
can provide a basis for sympathetic cooling of distant trapped clouds. 

PACS numbers: 32.80.Pj, 33.80.Ps, 42.50.Vk 



I. INTRODUCTION 

It is a well established fact, both theoretically and ex- 
perimentally, that light forces on atoms are substantially 
modified within resonant optical cavities P, S S 01 S IE 
0, 1^ • Possible experimental realizations range from sin- 
gle atoms or ions Jg, llOll in microscopic super cavities to 
several thousand |llUl2j| or up to a million atoms 0, 
in a high Q ring cavity. Applications of these systems 
include possible implementations of quantum informa- 
tion processing setups ITU , and controlled nonclassi- 
cal light sources jl^ as well as new possibilities for 
trapping and cooling of atoms and molecules. The basic 
physical mechanism in these setups can be traced back 
to the backaction of the atoms on the field. They act as 
a moving refractive index and absorber, modifying the 
intensity and phase of the intracavity field, which in turn 
governs their motion. This coupled dynamics is at the 
heart of cavity enhanced trapping and cooling. 

It is clear that if a single atom is able to change the 
field, it will influence other atoms in the same field ir- 
respective of their distance. This introduces long range 
atom-atom interactions, which are widely tailorable by 
suitable choices of cavity geometries and operating con- 
ditions. On one hand these interactions are useful and 
can be used to implement bipartite quantum gates |l4j . 
On the other hand they play a decisive role in the scaling 
properties of cavity enhanced coohng 0, 0| . For per- 
fectly correlated atoms, the change of refractive index 
induced by one atom can be compensated by a second 
atom, so that the effective atom-field back reaction can 
be strongly reduced. For several atoms in a ring cavity 
this effect only allows a weak damping of relative mo- 
tion, while the center of mass motion is strongly damped 
[20| . This model is closely related to the so-called CARL 
laser, where the kinetic energy of an atomic beam leads 
to gain into the counterpropagating mode of a single side 
pumped ring resonator [21|. New effects were also found 
in the study of the coupling of two Bose-condensates in 



a cavity |22j |. 

Several limiting cases for N atoms commonly inter- 
acting with a cavity mode have already been studied. 
For the case of N strongly trapped atoms in a standing 
wave cavity mode, it is possible to derive a set of cou- 
pled equations for the total kinetic and potential energy 
as well as the field amplitude |2^, which exhibit col- 
lective, damped oscillations ending in highly correlated 
steady states. This approach, however, does not give 
much insight into the details of the individual dynamics 
and correlations. In the opposite limit of N untrapped 
atoms moving in the cavity field, numerical simulations 
show little infiuence of atom-atom correlations and cool- 
ing proceeds independent of the atom number for 
proper rescaling of the cavity parameters. 

Recently an approach for several atoms in a single 
mode cavity has been developed, which concentrates on 
the effect of the — 1 other particles on the cooling 
properties of a single one 19]. This, in principle, makes 
it possible to study the combined optical potential and 
friction forces. It has been recently proposed theoreti- 
cally and confirmed experimentally |25, that if the 
atoms are pumped directly from the side (as opposed to 
pumping the cavity), the buildup of spatial correlations 
within a cloud of trapped atoms can lead to superradi- 
ant light scattering and enhanced cooling behavior. The 
theoretical results in this case are based on numerical 
simulations of the semiclassical equations of motion for 
a large number A^ ^ 1 of particles pSl- This clearly 
demonstrates collective effects, but does not give much 
quantitative insight in the buildup and role of atom-atom 
correlations. 

The central goal of the present work is to study the 
basic physical mechanisms responsible for the motional 
correlations and to develop quantitative measures of the 
established steady state correlation. For this we restrict 
ourselves to the simplest nontrivial example, namely two 
atoms strongly coupled to a single standing wave field 
of a cavity. The energy loss is compensated by cavity 
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tion picture read 0: 



FIG. 1: The setup. The motion of two two-level atoms in a 
far-detuned, high-Q optical cavity is studied. The cavity is 
pumped by a strong laser beam, almost on resonance with 
the cavity mode. Atomic motion is discussed along the cavity 
axis only. 



pumping, and large detuning from the atomic transition 
is taken to ensure low atomic saturation. Moreover, the 
motion of the atoms is only followed along the cavity 
axis. As we will see, this contains most of the essential 
physics but still allows us to derive analytical expressions 
for many relevant quantities. Most of the analytical re- 
sults are valid for the more general A^-atom case. 

The article is organized as follows: after presenting our 
model and the approximations used in section II, we ana- 
lytically discuss the central physical mechanisms present 
in section III. In section IV we quantify the results using 
numerical simulations, which are then analyzed in more 
detail in section V. 



II. THE MODEL 

Let us start by outlining the system which is shown in 
Fig. n We consider N — 2 two-level atoms with tran- 
sition frequency loa strongly coupled to a single mode 
of a high- finesse cavity with frequency ujc- The system 
is driven by a coherent laser field of frequency w and 
amplitude ry injected into the cavity through one of the 
mirrors. The model and its theoretical treatment follow 
closely Ref. 0. Coupling to the environment introduces 
a damping via two channels. First, the atoms sponta- 
neously emit with a rate of 27 into the vacuum outside 
the cavity. Second, the cavity photons decay with rate 2 k 
via the output coupler mirror of the cavity. The atoms 
can move freely in the cavity, however, for the sake of 
simplicity, their motion is restricted to the cavity axis 
(dashed hne in Fig.QJ. 

Applying the standard Born-Markov approximation 
the dynamics are governed by a quantum master equa- 
tion. 
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The atomic and cavity field detunings are defined as 
— Lo — to A and Ac = uj — wc, respectively. The an- 
nihilation and creation operators of the cavity field are 
d and d\ while and aj. are the lowering and raising 
operators of the A:-th atom. The Hamiltonian consists 
of the motional and internal energy of the atoms, the 
self-energy of the cavity field, the classical (laser) pump, 
and the Jaynes-Cummings-type interaction between the 
atoms and the field. The position dependence of the cou- 
pling constant is due to the spatial structure of the field 
mode: g{xk) = 50 /(x^), where in our standing- wave cav- 
ity /(xfc) = cos{kcXk), with kc = 2'k/X being the cav- 
ity wavenumber. The Liouvillean operator includes the 
effect of cavity losses and of the spontaneous emissions 
on the combined atom-cavity field density operator g. 
This latter is given by the last term, where the integral 
goes over the directions of photons spontaneously emit- 
ted by the atomic dipole, having expected wavenumber 
kA = i-^a/c and angular distribution A'^(u). 

We consider cold atoms but with a temperature well 
above the recoil limit ksTrec = ^k\/{2M), where M is 
the mass of one atom. In this limit the atomic coherence 
length is smaller than the optical wavelength and the 
position and momentum of the atoms can be replaced by 
their expectation values and treated as classical variables. 
We still keep the quantum nature of the internal variables 
d and (ifc. Moreover, if the atoms move much less than a 
wavelength during the equilibration time of the internal 
variables 



V < Ak, A7, 



(3) 



we can adiabatically separate the "fast" internal from 
the " slow" external dynamics as in standard laser cooling 
models k3. 



p^-'-[H,p\+Cp. (1) 

Using rotating-wave and dipole approximations the 
Hamiltonian and the Liouville operators in an interac- 



A. The internal dynamics 

For given positions of the atoms the internal 
atomic dynamics can be rewritten in the form of quan- 
tum Langevin equations. For low saturation, i.e. when 
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((7^(t) ^ 1, we can approximate the operator by —1/2 
(this is cahed bosonization of the atomic operators). The 
resulting Heisenberg-Langevin equations then reduce to 
the foUowing set of coupled linear differential equations: 



—a = (iAc - K)a + ^ g*{xk)ak + V + ^ (4a) 

k 



— (7fe = {iAA - 7)0"*: - 5(xfc)a + Cfc- 
at 



(4b) 



The noise operators ^ and Cfe appear as a result of the cou- 
pling to the external vacuum through the cavity mirrors 
and through spontaneous emission. They contain the an- 
nihilation operators of the external vacuum modes, and 
therefore give when acting on the environment's state. 
Their second-order correlation functions are as follows: 

{i{ti)i+it2))=2KS{h-t2), (5a) 

{Uti)C+it2)) = 2jSk„Ati - t2), (5b) 

while all other correlations vanish. 

The steady state expectation values of the inter- 
nal variables a and ak obtained from the Heisenberg- 
Langevin equations Q then read: 

7 - 
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Here V is the reduced determinant of the Bloch matrix, 
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Since the factor 1/T>' appears in both expectation values 
and in later formulae as well, it is worthwhile to rewrite 
it to reveal its resonance structure with respect to the 
cavity detuning: 

(8) 



v (zAa-7) i(Ac-c/)-(« + r)' 



where 



and 



U = 



Ai+7^ 



C/o5]/2(xO, (9) 
I 



and we used g{x) — gof{x). It is clearly seen that each 
atom broadens the resonance at most by Tq and displaces 
it by Uq. 



B. The external dynamics 

The motion of the atoms is governed by the force op- 
erator. 



''k = ■^[P'-^] = *^ V.g(xfe) o-]^^a - Vg*(xfc) a^o-fe 



Since is normally ordered, its expectation value is eas- 
ily obtained upon insertion of the stationary solution © 
of the internal variables. For a moving atom this expres- 
sion is only approximately valid: there will be a time lag 
in the internal dynamics with respect to the atom's cur- 
rent position, and hence we will include corrections to Ffe 
to first order in the atomic velocities. 

The slow evolution of the centers of mass of the atoms, 
smoothed out on the timescale r w max{ 1/k, I/7} is 
described by the coupled Langevin equations: 

ifc = Pk/M, (12a) 



N 



(12b) 



In these equations fk = i^'^kj are the vectors giving the 

steady state v = contribution of the force, while (3^^ 
are the tensors describing the first order corrections to 
the force acting on atom k. Hfc denotes the Langevin 
noise forces due to photon recoil. They correspond to 
random kicks along the cavity axis with zero average and 
second moments given by (3^3^) = 2?/cm- Note that the 
(6a) matrix Vkm, representing the strength of the Langevin 
noise, depends on the time-varying atomic positions. It 
(6b) represents the quantum fluctuations of the force due to 
the fact (^Fk o / (^Fk'^o {^F,„y It is by the addition 

of the noise terms that we tailor our classical force 
to give the same second-order expectation values as its 
(7) quantum comiterpartQ. 



III. INTERACTION CHANNELS 

The Hamiltonian |(2Jl contains no direct coupling be- 
tween the two atoms: these only arise indirectly due to 
coupling to the same field mode. Interestingly the atom- 
atom interaction appears in all the three types of forces 
present in classical equations of motion H12() . First, the 
steady state force ffc depends on the positions of both 
atoms via the steady state intensity. Second, not only 
does the friction coefficient on one atom depend on the 
position of the other, but the friction matrix has off- 
diagonal terms as well. This means that apart from or- 
dinary viscous friction (pfe oc Vfc) a strange phenomenon, 
which we call cross- friction (p^ oc v;, for Z ^ k) is also 
present. Here the velocity of one atom influences the fric- 
tion experienced by the other atom. Third, the Langevin 
noise term on one atom has an expected magnitude influ- 
enced by the position of the other, and the noise terms 
Hfe are directly correlated as well. Hence we get joint 
"kicks" on both atoms leading to correlated motion. In 
the following we will analyze these interaction channels 
(11) in more detail. 
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A. Steady state force f 

Formally the expectation value of the force opera- 
tor looks very similar to the case of free-space Doppler 
cooling: 

h = -h--^^{a^d)Vu9\^k). (13) 

However it additionally depends on the positions of the 
other atoms via the cavity field intensity. This dipole 
force is conservative and can be derived from a potential. 
The potential looks more complicated than in free space 
as it contains the dynamical nature of the cavity field, 
but still can be given in closed form 0| : 

F^/^^'f atan^--^/^^+,^'^^^-'). (14) 

To show the effects of dynamic field adjustment, we 
plot this potential for two typical cases in Fig.jS] For the 
experimental parameters used in an experiment at MPQ 
in Garching [23| (Garching parameters), k = 7/2, go = 
57, and for a detuning A^ — — 5O7 (upper graph) the 
effective interaction between the atoms is relatively weak 
and the potential resembles the familiar "egg-carton" sur- 
face proportional to sin^(fccxi)-t-sin^(fcca;2). For a some- 
what stronger atom-field coupling = 2O7 (lower graph 
in Fig. 121) the atomic interaction is quite obvious and the 
shape of the potential of the second atom strongly de- 
pends on the position of the first atom and vice versa. Ba- 
sically in this second case either both atoms are trapped, 
or both are free. 

The peculiar go-dependence of the interaction and the 
trapping effects can be understood physically by look- 
ing at formula The atoms see each other through 
the cavity field (a^a). As we saw in formulae (|8I9I10|I . 
the field is in resonance when U (xi , X2 ) is approximately 
Ac ± (r + k) . Each atom can shift U in this far-detuned 
case by approximately Uq ~ g^/ IS.a, whereas the cavity 
linewidth is approximately K + -fglllS,\. Using the MPQ 
parameters < k, and therefore the back-action of the 
atoms on the cavity field is weak. In the large-go case, 
Uq > K, but Fq < At, meaning that the atoms can shift 
the cavity resonance significantly more than a linewidth. 
Hence if one of the atoms leaves its trap it will shift the 
cavity out of resonance and cause the other atom to be 
released as well. Moreover, the amplitude of the force 
is proportional to the cavity field, which in the large-go 
case decreases faster with the distance from the trapping 
point. This implies that less work needs to be done to 
free an atom: the potential is smaller, despite the same 
maximum light shift. 

Let us now analyze the interaction in a more quanti- 
tative way. We start from the limiting case of neglecting 
the atomic back-action on the field (i.e. we assume a con- 
stant field intensity). This corresponds to the force Ijl^fl 




FIG. 2: Potential as a function of atomic positions The po- 
tential 11411 is plotted as a function of positions of the two 
atoms. In a) the Garching parameters (k = 7/2, go = 67) 
are used, in b) go was increased fourfold. In the first case, 
the potential is well approximated by a sum of two single- 
particle potentials. In the second case, it is seen that either 
both atoms are trapped or both are free. The trap is deeper 
for smaller couplings. 



one obtains in the limit of very large atom-field detuning. 
We can now find corrections from the cavity-mediated 
interaction to this force by expanding the potential H14() 
in a power series. Setting the driving field to resonance, 
A(7 — Uq — K, the potential Eq. (|14|l to second order in 
7/A^ reads 




To first order in the small parameter we have a sum 
of single-atom potentials, giving the "egg-carton" shape. 
The corrections to this are given by terms of higher or- 
der in our expansion parameter, of which we give the first 
nontrivial term here. Note that since it is not simply the 
distance of the atoms upon which the potential depends. 
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the interatomic force between them - as can be read out 
from the above formula - is not a "force" in the sense of 
Newton's third law. 



B. Forces linear in velocity: friction and 
cross- friction 



To lowest order in the adiabatic separation of the in- 
ternal and external dynamics we used the steady state 
values of the internal variables for fixed positions of the 
atoms to calculate the above potentials. As a next step 
we can include corrections for a and the ak linear in the 
velocity of each atom, which should be valid for low 
velocities. As described in I'] this leads to a friction ma- 
trix /3fc„ as first order correction to f^. 

We obtain the following explicit formula for the friction 
matrix, 



/3/cm 2Wfcg(xfc) o V™g(x„ 



+ Wfc52(Xfc)oV™52(x,n 



2\V'\^ \V'^ 



2(1 + X) 



(^(*A^-7)^-^.g(x,f j 



+ {1 + 3X)V' 



(17) 



where x = {i^A + l)/{iAA — 7) is a complex factor of 
unit modulus, which for large atomic detuning, ^ 7, 
becomes approximately x ^ 1- This formula differs 
somewhat from that obtained by Fischer et al. by 
a slightly different approach. The most important differ- 
ence is that we find a matrix symmetric in the indices 
km. This is important because it is these off-diagonal 
terms that couple the velocities of the atoms, and have 
a decisive influence on the buildup of correlated motion. 
We defer detailed discussion of the formula (|17|l to Sec- 
tion |V| 



C. Random forces due to quantum noise 

Spontaneous emission and cavity decay introduce 
quantum noise into the atomic motion. These heat up 
the system and generally tend to decrease motional cor- 
relations of the atoms. Following the line of reasoning 
briefly mentioned at the end of section ^ and discussed 
in more detail in Q , we can calculate the influence of the 
noise operators ^ and C, of eq. Q on the dynamics. For 
N atoms we arrive at the following simple formula: 



■n 



Vkni = 2^1 Vfeg(xfe) o V™5(x,„)-;^^74m 



2h^Vk9{-^k) o V„ig(x„ 



2 kA^ + J Ac 



\V'\ \V'\ 



This is a simple extension of the corresponding formula 
for one atom given in Let us remark here, that the 
diagonal part of this diffusion matrix Vkm has been also 
found by Fischer et al. T^. Surprisingly one also obtains 
off-diagonal terms, which have not been considered be- 
fore. These terms lead to correlated kicks on the atoms, 
which can .add to the atom- atom correlations, rather 
than destroying them. 

Spontaneous emission adds recoil noise, which gives an 
extra term to the noise correlation matrix of the form: 



km 



-1- 



(19) 



(18) 



Here kA — ^a/c is the expected wavenumber of the 
emitted photons and is the correction factor coming 
from the spatial distribution of the photons, in our case 
V? — 2/5. As expected, spontaneous emission, being a 
single-atom process, induces no correlations between the 
atoms. 



IV. NUMERICAL SIMULATIONS OF THE 
CORRELATED ATOMIC MOTION 

Having discussed the qualitative nature of the com- 
bined atom field dynamics, we now turn to numerical 
simulations for some quantitative answers. We numeri- 
cally integrate the Langevin equations (|12|) . varying the 
ratio of the cavity loss rate to the linewidth of the atom 
k/7 as well as the relative coupling strength go/7- Note 
that the relative magnitude of radiation pressure and 
dipole force can be changed by varying the detuning be- 
tween pump frequency and the atomic resonance. As we 
are interested in cooling and trapping the atoms we fix 
the cavity frequency at Ac = NUq — k to ensure effi- 
cient cooling ll|,|28|. The pump power is always chosen 
to keep the atomic saturation low and approximately at 
(((T^o-fe) < 0.1) at all times. 

Typical trajectories of atom pairs are shown in Figure 
13 The atomic positions (xi (i„), a;2(tn)) are plotted at 
regular time intervals i„ . For better visibility the coordi- 
nate is always measured from the nearest trapping point. 
In the left column the detuning was chosen (A^ = — 5O7), 
with the cavity parameters of the MPQ group at Garch- 
ing |23| (k = 7/2, go — 57). The first 50/^s are displayed 
in the figure on the top, and the first 3ms in the one on 
the bottom. Both atoms localize during the first 50/j,s 
to within 1/4 of a wavelength around respective trap- 
ping centers, a sign of cooling and trapping by the cav- 
ity. In the right column, the detuning is chosen much 
larger (A^ = — 10''7) with stronger atom-field coupling 
(k = 7/10, go = IO7). In this case the relative impor- 
tance of spontaneous emission is strongly reduced. The 
cooling in the cavity is faster and both atoms reach a 
steady state rapidly. 

The appearance of a circular structure is the striking 
feature of the right column. This indicates that the mo- 
tion of the atoms is correlated. Both atoms move si- 
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FIG. 3: Typical trajectories in coordinate space. In the left 
column, in a) and b), the parameters of the Garching group 
are used, in the right column, in c) and d), we considered 
a better resonator (k = 0. 17,30 = IO7) and larger detuning 
(Aa = — IOOOO7). In both cases the first 50/.is, in a) and c), 
and the first 3ms, in b) and d), are shown. The coordinate is 
the distance from the nearest trappping point. 




t[ms] 



10 



t[msl 



012345 56789 1 

t[ms] t[ms] 



FIG. 4: The phases of the atomic oscillations and the relative 
phase. The time evolution of the oscillator phases of the first 
(upper row) and second (middle row) atoms is shown, with 
the phase difference (bottom row). On the left the parameter 
set of the Garching experiments is used, on the right a better 
cavity is taken with larger atomic detuning, as in the text. 
The phases of the atoms evolve too fast on this timescale, 
only noise is seen. The phase difference is slower, with the 
"improved" parameter set it appears to stabilize at -f 90° and 
-90°. 



nusoidally about their respective trapping points, with 
some noise but in such a way that the relative phase of 
the two oscillations is likely to be +90° or —90°. 

To quantify the correlation between the atomic oscil- 
lators we define their oscillator phases. This is computed 
in the simulation using the trap frequency, which is: 

entrap = ^2n|A^|(ata) kl/M, (20) 

if both atoms are well trapped. Here kc is the resonator 
mode wavenumber, and (o'^o') is the saturation of either 
atom at the trapping point. This formula can be derived 
by expanding the potential H14|) . and substituting our 
particular choice of r/ and Ap. 

The measured time evolution of the oscillator phases 
for the MPQ parameters (left column) and the "im- 
proved" parameters (right column) are shown in Fig. 0] 
For each parameter set, representative runs shown in 
Fig. 13 are used, and the oscillator phases of atom 1 (up- 
per row), atom 2 (middle row) and the phase difference 
(lower row) are displayed. The rapid oscillation of the 
atoms (the period is 2.9 /.ts for the Garching parame- 
ters and 0.17 /is for the idealized ones) means that the 
time evolution of the phases is too fast to be followed on 
the timescale shown, in the upper and middle rows only 
"noise" is seen. The phase difference, however, evolves 
more slowly. This effect is more pronounced for the sec- 
ond parameter set, where the dipole force dominates. 
Moreover, in this second case, the phase difference clearly 



stabilizes around -|-90° or —90°, with random jumps in 
between, as expected from the results shown in Fig. El 

What is left is to define a single number that quantifies 
the strength of the time averaged correlations. For this 
we sample the distributions of the oscillation phases and 
the phase difference over time to create the histograms 
shown in Fig. [S] As we expect, the distribution of the 
phases is relatively flat for both parameters sets. How- 
ever, there is a significant difference in the distribution of 
the relative phase: we find very pronounced peaks around 
-1-90° and —90° for the second parameter set. 

The asymmetry in these peaks is a numerical artefact 
due to the finite sampling time. It is strongly diminished 
if we average over several different initial conditions. A 
parameter that measures the magnitude of the corre- 
lation is the width W of these peaks, defined through 
= (|A0| — 90°)^, where the overbar denotes averag- 
ing over time and over several trajectories with differ- 
ent initial conditions. Subtracting this width from the 
width of the flat distribution we get the signal strength 
S = 51.96° — W , which can be obtained directly during 
the simulation. 

Let us point out here, that this measure S of motional 
correlation is useful only if the atoms are well trapped. In 
fact, fast atoms freely moving along the lattice generate 
peaks in the single atom phase distributions of (f>i or (/)2 
around 0° or 180°. We therefore monitor the single atom 
distributions simultaneously and measure the widths of 
the peaks around 0° and 180° in the same way as we 
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FIG. 5: The distribution of the phases (f&i, $2^ and the rela- 
tive phase (^12). For the Garching parameter set, in a), the 
distributions are flat. Slight peaks in $1 and $2 at 0° and 
±180° signal fast and free atoms. With the improved param- 
eters, in b), the trapping is better. Strong peaks at ±90° in 
$12 indicate correlation of atomic motion. The strength of 
these peaks is quantified by their width, both for the phase 
(Noise) and the phase difference (Signal). 
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FIG. 6: Correlation vs k and go at various detunings. At 
small detunings |Aa| < IOO7 the atoms are cooled but not 
trapped, correlation cannot be measured. At higher detun- 
ings trapping is good enough, and the correlation strength 
measured as described in the article is shown in with shades 
of gray. 



did for the signal S. If the "Noise strength" is too large, 
(> 10°, an arbitrary value), we declare our correlation 
measure unusable. 



A. Scanning the parameter space 

Having defined a suitable measure of correlation be- 
tween the motion of two atoms in the same cavity, we 
are in the position to quantitatively investigate the pa- 
rameter dependence of this phenomenon. To this end, 
we ran the simulation program for cavity decay rates of 
1/50 7 < K < 57 and couphng strengths 7 < 50 < IOO7, 
at atomic detunings of —10^7 < A^i < — 5O7 (at de- 
tunings of higher magnitude the atoms move too rapidly 
and adiabaticity does not hold). The results are plotted 
in Fig. El 

At atomic detunings of IOO7 or less, the atoms are 
cooled but generally not trapped by the cavity, except 
for a small region of go < 57 and k < 7. In that case the 
motion of two atoms will not be correlated. For detunings 
as large as 5OOO7, the cavity field traps and cools the 
atoms for any considered values of the parameters. The 
correlation becomes apparent for go > 40k, and grows 
weaker if go is further increased. 



B. Correlation and cooling 

A decisive advantage of cavity cooling of an atom is the 
fact that it needs little spontaneous emission '2^ . It has 



been argued that this is a single atom effect and correla- 
tions established between the atoms' motion will decrease 
efficiency or even turn off cooling for a trapped ther- 
mal ensemble. On the other hand, in some simulations 
for very weakly bound atoms, this effect seemed not to 
play any role [l8| . 

In our model we can now study the effect of correla- 
tions on cavity cooling in a very controlled way for a large 
range of parameters. We performed several runs of our 
simulation with different cavity parameters and at dif- 
ferent detunings, comparing the equilibrium temperature 
for one and two atoms in the cavity. Our results are plot- 
ted in Figs. 13 (one atom) and|Hl(two atoms). All shades 
of gray show temperatures below Doppler temperature 
hrf, while black denotes temperatures above that value. 

Putting one atom into the cavity (Fig. [TJ, we find 
that sub-Doppler cooling is achieved in the good-coupling 
regime go > ^,7. In that regime, the final temperature 
does not depend on go, but is approximately proportional 
to K, as expected from previous work |26l l29j|. Here one 
has to be cautious of the results, since if the atom's ki- 
netic energy is only a few times the ground-state energy 
of the harmonic trap potential, the validity of the semi- 
classical approximation can be questioned. 

Interestingly, the temperature plots look different if we 
load two atoms into the cavity (Fig. At moderately 
large detunings, Aa — — 5O7 and — — IOO7, temper- 
atures are the same as in the single-atom case. In these 
cases the atoms are cooled, but not trapped, by the cavity 
field. At larger detunings, = — 5OO7 and larger, we 
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FIG. 7: One atom in the cavity: final temperature vs. k and 
go, for different atomic detunings. Temperature is shown in 
units of h'y with shades of gray. Black denotes temperature 
above the Doppler limit of ^7. 
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FIG. 8: Two atoms in the cavity: final temperature vs. k and 
go, for different atomic detunings. Temperature is shown in 
units of h'y with shades of gray, as in Figure |7| For low atomic 
detuning |Aa| < IOO7 the temperature is the same as in the 
one-atom case. For large atomic detuning new structures ap- 
pear in the plots. 



see regions in the plot where "extra heating" compared 
to the single- atom case is observed. This effect is most 
prominent at extremely large detunings — — 5OOO7 
and A A = — 10'*7, where for k < 7 a new structure ap- 
pears in the temperature plots. This points to the highly 
reduced efficiency of cavity cooling: whereas for a single 
atom a better resonator (lower k) implies lower final tem- 
perature, if there are two atoms in the cavity, decreasing 
K can increase the temperature. 

Comparison of the temperature plots in Figs. [7| and |S1 
with the plots of the correlation strength in Fig.Elreveals 
strong similarities. Indeed one sees that the excess heat- 
ing caused by the presence of the other atom coincides 
with the buildup of correlations in the motion. In other 
words, the correlation established between two atoms re- 
sults in a loss of efficiency of cavity cooling, with final 
temperatures pushed up to the Doppler limit. 



C. The origin of correlation: numerical tests 

As we have seen above, the motion of the atoms be- 
comes correlated due to the cavity-mediated cross-talk. 
This interaction occurs via the force, via the friction, and 
via the diffusion as well. We would like to know which of 
these interaction channels is responsible for the correla- 
tion. In our simulation we can conveniently answer this 
question if we artificially weaken only particular channels 
of interaction and measure the correlation. 

The interaction can be eliminated from the determinis- 
tic force using the approximation (|15|l . Friction and dif- 
fusion contain interaction at various levels. On the one 



hand, cross- friction and cross-diffusion are direct interac- 
tions between the atoms. On the other hand - through 
the determinant V - the positions of both atoms influ- 
ence the friction and diffusion constants for the other 
atoms. This parametric interaction was not analyzed, 
but does not seem to play a prominent role. Both cross- 
friction and cross-diffusion can be eliminated by simply 
suppressing the off-diagonal terms in the friction and dif- 
fusion matrices. 

The elimination of the interaction can be made contin- 
uous with mixing parameters < yF,yp,yv < 1, giving 
the force, the friction, and the diffusion with the follow- 
ing formulae: 

F' ^yFF+{l-yF)Fo, (21) 

P'km = Vf3Pkr,i + 4m (1 " yp)Pkrn; (22) 
2^fem = yv'Dkrn + 4m(l " yv)'Dkm- (23) 

We show an example of what this gives for detuning 
= — 5OOO7, decay rate k = 0.67 and coupfing 
go = 3O7 in Fig. |5| The results of these investigations 
are quite unanimous. It can be seen that linearization 
of the potential has no systematic effect on correlation 
strength, and correlations are slightly enhanced if the 
noise is decorrelated. If cross-friction is eliminated, how- 
ever, all correlations disappear. We can therefore con- 
clude that the dominant effect leading to correlated mo- 
tion is related to cross-friction. 
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FIG. 9: Correlation for various interaction channel strengths. 
The interaction via the conservative, the friction, and the dif- 
fusion forces is weakened by the mixing ratios yp, yfj and yo, 
as described in the text. The correlation strength was then 
measured using our correlation parameter. Cross-friction is 
seen to play the dominant role in the establishment of corre- 
lation. 



V. THE ORIGIN OF CORRELATION: 
ANALYTICAL EXPANSION FOR DEEP 
TRAPPING 

The emergence of motional correlations between 
trapped atoms can be examined analytically using the 
formulas presented in Section Hill Both friction and dif- 
fusion matrices have an "isotropic" term proportional to 
Skm , which acts on each atom separately, and an "interac- 
tion" term, containing cross-friction and cross-diffusion. 
The latter can be rewritten in the following form, using 
the notation gk — g{^k)- 



-^PiM) 




(24) 



where the prefactor P depends on the positions of all 
atoms symmetrically. 

For well-trapped atoms a;^ <C A (xk denoting the dis- 
tance of atom k from the nearest trapping site), the cou- 
pling constants can be approximated as 



gk = gixk) = go cos{kcXk) « 50 l - 7: fc 



1 



2 2 



(25) 



Substituting these into the interaction parts of the ma- 
trices (OH yields: 



P{{^i})gtkt,r^ 



(26) 



where the "radius" , r — xf) ^ , is the distance of the 
system from the origin in the coordinate space. This ma- 
trix is a projector to the vector (xi, . . . ,xjv), and hence 
the interaction terms affect only the "radial" motion. 



In the well-trapped case we can write the friction and 
diffusion matrices to second order in the coordinates xi 
as: 



Pkm — Po^km + PlXk^m 
T^km = T^oSkm + ViXkXr, 



The following notation is used: 



r 2,2 2Aa 



l3o = 2h oignkr^Ty ^knXu 



2 

2?o = 2?i ^IQokn 

IP' I' ° ^ 



+ 7^ 
2;.2 t 2 _^ 3^2 ^2 



2?o(l-3x) + 



2(l-x)((*A^-7)'-E3(^')')]| 
V,^2n^^gtkhi^^^^^, 



(27) 
(28) 



(29) 
(30) 

(31) 
(32) 
(33) 



where V'q stands for the Bloch determinant V defined in 
eq. (O, evaluated at the trapping site. Near the trapping 
site all Xk go to zero, and only remains finite. This 
prevents the atoms from stopping completely: they are 
"heated out" from the trapping sites themselves. At some 
distance from the trapping sites the interaction terms 
become important as well, in the coordinate space this 
induces extra diffusion and friction in the radial direction. 
Provided go is high enough, radial friction is enhanced 
much more than radial diffusion, leading to a "freezing 
out" of the radial mode, i.e. a motion at approximately 
constant distance from the origin. 



VI. CONCLUSIONS 

As pointed out in several previous papers, the motion 
of particles in a cavity is coupled quite generally through 
the field mode. This implies the buildup of motional cor- 
relations, which we have investigated here in more detail. 
In general we found that in steady state these correla- 
tions are hard to see directly and it is difficult to find 
good qualitative measures to characterize them. Mostly 
they are strongly perturbed by various diffusion mech- 
anisms. However, they still play an important role in 
the dynamics and thus can be observed indirectly. As 
one consequence they can lead to a significant change 
(increase) in the steady temperature, which directly re- 
lates to trapping times and localization properties. This 
poses limits to cavity induced cooling for large particle 
numbers. 

From the various mechanisms at work, cross-friction 
turns out to be the most important. It creates corre- 
lation and leads to a fast thermalization of two distant 
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particles without direct interaction. Hence, this mecha- 
nism should prove vital for the implementation of sym- 
pathetic cooling of distant ensembles coupled by a far 
off-resonant cavity field. In order to get efficient cou- 
pling the two species should have comparable oscillation 
frequencies, so that correlation buildup and thermaliza- 
tion is fast. The finesse of the cavity should also be large 
(long photon lifetime). 

The second important coupling mechanism, which 
works via the joint steady state potentials, was suggested 
for use in the implementation of conditional phase shifts 
[T^ . It can be viewed as a cavity-enhanced dipole-dipole 
coupling. Although this contribution will become more 
important for larger atom-field detunings, the conditions 
for this part to dominate the dissipative cross-friction 
seem rather hard to achieve in practice. Finally we found 
that the noise forces acting on different atoms contain 



nonlocal correlations too. These are particularly impor- 
tant for large detunings and relatively low photon num- 
bers, where spontaneous emission is strongly reduced. As 
a result they could seriously perturb bipartite quantum 
gates in cavities. 
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